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Abstract. The problem of an electromagnetic wave scattered from a random 
medium layer with rough boundaries is formulated using integral equations which 
involve two kinds of Green functions. The first one describes the wave scattered 
by the random medium and the rough boundaries, and the second one which 
corresponds to the unperturbed Green functions describes the scattering by an 
homogeneous layer with the rough boundaries. As these equations are formally 
similar to classical equations used in scattering theory by an infinite random 
medium, we will be able to apply standard procedures to calculate the coherent 
field. We will use the coherent potential approximation where the correlations 
between the particles will be taken into account under the quasi-crystalline 
approximation. 



§ asoubret@hms.harvard.edu 



Electromagnetic wave scattering from a random layer with rough interfaces I 



2 



1. Introduction 

Many studies on electromagnetic waves scattered by a random medium layer with 
rough boundaries have been reported in recent years PElElEElEllZllHlini EH 
ITT1 IT21 IT~51 IT31 ITS] . Rigorous numerical methods have been developed |21 E| but 
are computationally intensive or limited to 2D geometry. Most often, the radiative 
transfer theory is used for the volumetric scattering with the Kirchoff or small- 
perturbation method for imposing the boundary conditions 0] |5J [7| |SJ H E| • 
This method is well suited to compute the scattered intensity but is based on 
phenomenological considerations. Thus, analytical theory has been developed in order 
to describe the coupling between the random medium and the rough boundaries. 
Furutsu |13l 114) formulates the rough surface scattering problem with Dyson and 
Bethe-Salpeter equations which permit treating the random medium and the rough 
boundaries on the same footing. Unfortunately, this approach is formal, and the 
relationship between the radiative transfer theory and the classical rough surface 
scattering theories [T7| HH1 III E01 HH is not straightforward. Mudaliar [TUIITTI IT2] 
uses integral equations where the rough boundaries are treated under a perturbative 
development. He shows that the intensity verifies a "generalized" transport equation. 
If this approach is more numerically tractable than Furutsu's, the expressions obtained 
are still involved. This is due to the choice of perturbative development to describe 
the scattering by the rough surfaces. In this paper, we show that we can obtain the 
general expression, whatever the choice of the scattering theory used at the boundaries, 
in introducing the scattering operators of the rough surfaces [211 . Furthermore, in 
separating the surface and the volume scattering contributions with the help of Green 
functions, we will be able to use well developed analytical theories of waves scattered 
by an infinite random medium [23 ESI 113 12H I2S1 EZl EE] In this paper, 
which is the first part of a series of three papers, we investigate the coherent field 
scattered by the rough surfaces and the random medium. The contribution of the 
random medium will be taken into account in introducing an effective permittivity, 
which is calculated under the Quasi-Crystalline Coherent Potential Approximation 
(QC-CPA) |lfil I17j . The contribution due to the rough surface will be given by the 
average of the scattering operators [2]. The calculation of the incoherent fields will 
be the subject of the following papers, where the derivation of the radiative transfer 
equation will be detailed, and the particular case of strongly diffusing random media 
will be treated using a vectorial diffusion approximation for Rayleigh scatterers. 

2. Geometry of the problem and formulation 

The geometry of the problem is shown in Figure ^ Volumes Vq and V% are 
homogeneous media with permittivity 6q(uj) and 62(0;). For simplicity we suppose 
that eo^) is a real positive number. The random medium V\ is made of spherical 
scatterers of permittivity e s (u) in a background medium of permittivity ei(w). The 
boundaries are described by the random functions z = hi(x) and z = —H + Ii2(x). 

In the following, we consider harmonic waves with e _1 u 4 dependence. For a 
point source located at r$ = Xq + zq e z in the medium Vo, the field scattered by the 
rough surfaces and the random medium at the point r = x + z e z in the media Vq , 
Vi, V2 are, respectively, given by the dyadic Green functions 16, 29 G sv (r , r , u>), 
G sv (r,r a ,u>), G sv (r,r ,uj). Here, for Ggy (r ,r ,u), the upperscripts a,a are, 
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Random medium 



Figure 1. Random medium with rough boundaries. 



respectively, the receiver location and the source location. These Green functions 
satisfy 



• Propagation equations : 



■=00 



■00 



V x V x G sv (r,r Q ,uj) ~e (uj)K vac G sv (r,r 



S(r-r )I, 



=10 



(1) 
(2) 

V x V x Gs V (r 7 r ,u;) - e 2 (u;) K^ ac G 2 s V (r,r ,uj) = 0, (3) 

with K vac = uj/c vac the vacuum wave number, and c vac the light speed in the 
vacuum. The permittivity ev(r,u>) inside the random medium V± is defined by : 



V x V x G sv (r,r ,uj) - e v (r,uj)K; ac G sv {r,ro,u) = 0, 

-20 , „ , , „o -^=20 



N 



e v (r,Lu) = ei(u) + > (e s (w) - £i(w)) 9 s (r - r,-) , 



(4) 



where rj., .. . , t*at are the center of the particles, and <d s describes the spherical 
particle shape : 

' 1 if IHI <r s 
if ||rl|>r a ' 



9.(r) 



(5) 



with r s the particle radius. 
• Boundary conditions on the upper rough surface: 

n sl ■ e (uj)G g V (r,ro,uj) = h sl ■ e^u) G sv (r,r ,u;), 

7=?°° i \ - 7=? 10 i \ 

n s i x G sv (r,r , lu) = n sl x G sv [r , r , w) , 



n s i 



V x G sv (r,r ,uj) 



n s i 



V x G sv (r,r ,u>)) 



n sl x V x G sv {r,r ,uj) = n sl x V x G sv {r,r ,uj, 



(6) 
(7) 
(8) 

(9) 
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where r = x + h\{x)e z , and h s \ is the exterior normal to the rough surface 
z = h\{x): 

e z - Vhi(x) 



nsl ~ (l + (V/ii(x)) 2 )V2 ' 
• Boundary conditions on the bottom rough surface : 

n s2 ■ ei(w) G sv (r,r , w) = h s2 ■ e 2 (w) G sv (r, r , w), 

7=710 , . . -=20 . . 

n s2 x G sv (r,r 0l lu) = n s2 x G sv (r , r , w) , 



w S 2 • 
n S 2 x 



=10 



V x G sy (r,r ,w) 

_ —10 . 

V x G sv (r,r ,uj) 



n s2 



,20 



V x G sy (r,r ,w)) 



n s2 x 



_ --=20 , . 

V x G sv (r,r ,w / 



(10) 

(11) 
(12) 

(13) 
(14) 



where r — x + [-H + h 2 (x)]e z , and n s2 is the exterior normal to the rough 
surface z = h 2 (x): 

-e z + Vh 2 (x) 



n s2 



(15) 



(1 + (V/12^)) 2 ) 1 / 2 ■ 

• Radiative conditions at infinity in the media V and V 2 . 

We also use Green functions where the source is situated in the medium V\. The 
fields in the medium Vq, V\, V 2 are given by the Green functions G sv (r,r ,u)), 

1 1 2 1 

G sv (r,r ,uj), G sv (r,r ,uj) which verify: 

• Propagation equations: 

V x V x G <> s 1 v (r,r )-e ^)K^ ac G O s 1 v (r,r o ) = 0, 

V x V x G^ v {r,r ) - e v {r,u) K^G^r,^) = S(r - r )I , 

V x V x Gf v {r, r ) - e 2 (w) if 2 ac S^yfo r ) = , 

• Boundary conditions on the upper rough surfaces: 



(16) 
(17) 
(18) 



h sl ■e 1 (w)G 1 s V (r,r ) = n sl ■ e (w) G°g V (r, r ) 
h sl x G S y(r 7 ro) = n sl x GsV( r , r o) , 



n 3 i 



n s i x 



V x G s \(r,r ) 
V x G"" y (r,r ) 



n s i 



V x G sv (r,r ) 



n sl x 



V x G sv (r,r ) 



• Boundary conditions on the lower rough surface: 

H 21 

n s2 • ei(w) G sv (r,r , lj) = n s2 ■ e 2 (uo) G sv (r, r , w), 

1 1 2 1 

n s2 x G sv (r , ro, to) = h s2 x G sv (r,r , u) , 

h s2 ■ V x G sv (r,r Q ,uj)) 



n s2 



n s2 x 



V x G sv (r,r a ,uj) 
V x G sv (r,r ,u) 



n s2 x 



V x G sv (r,r ,Lu) 



(19) 
(20) 
(21) 

(22) 

(23) 
(24) 
(25) 

(26) 



• Radiative conditions at infinity in the media and 2. 



In order to to separate the contribution from the rough surfaces and the random 
medium, we introduce the dyadic Green functions G s (r,ro,u>), G s (r,ro,u>), 
G s (r, ro, to), G s (r, r , u), G s (r, r , to), G s (r,r ,uj) which describe the scattering 
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by the layer with the rough boundaries but without the random medium. These 
functions verify similar propagation equations and boundary conditions as the Green 
functions G°gy (r, r , ui), where the permittivity ey{r,uS) due to the random medium 
is replaced by an effective permittivity e e (uj) in equations (|21 ll7(l and the permittivity 
€i(uj) is replaced by e e (cj) in equations HT1 H$l . This effective permittivity will 
be determined using the Coherent-Potential Approximation (CPA) with the Quasi- 
Crystalline Approximation (QCA) [H El EH E3 ESj ■ We will show in Section H 
how to write these Green functions with the help of scattering operators, which are 
common tools in scattering theory by rough surfaces |21| . 



3. Integral equations 

The previous system of differential equations with boundary conditions can be 
transformed into integral equations (161 1171 [2"§| . For a source in medium 0, we have 



7=700 

^sv 


7=700 




■ v 11 


—10 

' ' 




■ v 11 


7=711 


•V 11 


7=710 


(27) 


7=710 

^sv 


7=710 




■ v 11 


—10 

■ ' 




■ v 11 


7=711 




—10 

■ I 


(28) 


7=720 

^sv 


7=?20 
— Org - 




■ v 11 


7=710 




■ v 11 


7=711 


■ v 11 


7=710 


(29) 



and for a source in the medium 1 



with 



7=701 -*fll 

^sv — *~*S ~1 


01 W n -^11 
- Ixg • V ■ l* sv , 


(30) 




11 prll >r ll 


(31) 


r 21 -r 21 a 


- G s ■ V ■ G sv , 


(32) 


V (r,r ,uj) 


= S(r-r )V 1 (r), 


(33) 




= K vac(tv(r,uj) - e e (u))I, 


(34) 


g definition : 






[A-B}{r,r ) 


= / d 3 ri A(r,ri) ■B{r ll r ) . 

JVi 


(35) 



A direct demonstration of these equations involves integral theorems [33], but it is 
easier to invoke the uniqueness of the solution and verify a posteriori that the integral 
equations i|27l32f> satisfy the propagation equations and the boundary conditions. 
For example, to demonstrate that equation i|17fl is verified, we apply the operator 
V x V x — e e K 2 ac I on (|31|l . and using the propagation equation satisfied by G s with 
the definition in (|33I35(I . we obtain 

(V x V x -e e K 2 vaL l) ■ G S y(r,r a ) = (V x V x -e e K 2 vac l) ■ G s \r, r„) 



11 



l(VxVx-a„y)'G s V -Gsyfaro), 

= 5(r - r )I + V 11 ■ G'sy (r,r ), 

= S(r - r )I + K 2 ac (e v (r) - e c )5^(r,r ) , 

which is the propagation equation in (|17J) . By using the same procedure, we show that 
the propagation equations I|1I31 116118(1 and the boundary conditions 1(7191 112114(1 and 
<(2UI221 124126(1 on the rough surfaces are satisfied. The boundary conditions at infinity 
in media and 2 are specified by the choice of a retarded Green function for G s . 




However, due to the introduction of the effective medium e e , the boundary conditions 
<|Hl llll[T^ll23(l are not satisfied, and we obtain the following boundary conditions: 

n sl ■ e e (u) G sv (r,r ,u) = h sl ■ e (u) G sv (r,r ,v), (36) 
h s2 ■ e e (ui) G sv {r 7 r ,uj) = h s2 ■ e 2 (w) G sv (r ,r ,u), (37) 
h a i ■ e e (uj) G sv (r,r ,oj) = h sl ■ e (w) G sv (r,r ,uj), (38) 

11 21 

h s2 ■ e e (uj) G sv (r,r ,uj) = n s2 ■ e 2 (u) G sv (r ,r ,uj). (39) 

We see that in the left-hand side of equations H36I39|I . the permittivity is not ei(u), 
as it must be, but is e e (uj). If we had defined the Green function G s ° describing the 
scattering by a homogeneous medium (with rough boundaries) with the permittivity 
ei, the problem would not exist. But if we want to use the Coherent-Potential 
Approximation, we must introduce this effective permittivity. We might go around 
this problem in changing the definition of the Green function G^ where a small 
layer of arbitrary small thickness h with permittivity t\ is added along the rough 
boundaries. (See Figure @.) These Green functions verify the boundary conditions 
(HO El 013 l23l) where the permittivity e\ is not replaced by e e due to the added layers 
along the boundaries. Therefore, with this definition, equations (|27I32I) verify the 
boundary conditions El IE3 ■ However, the propagation equations (0 H7I) are 
not satisfied since the added layers produced new contributions. But as we can choose 
the layers' thickness as thin as we want, we can neglect the effect of these layers on 
the propagation equations. In the following, we won't take care of these boundary 
condition problems, and we will suppose that equations (|27I32|) are solutions of our 
problem. The integral equations I|27I32|I are the key point of our approach. We see 
that to calculate the field in medium or 2 (when the source is in medium 0) with 

equations l|27|l and (|29|l . we first need to determine the Green function G sv , where 
the source and the receiver are in medium 1. This can be done with equation (|3f li 

where the only unknown is G sv . If the permittivities of the medium 0, 1, and 2 and 
the effective permittivity were equal (e e = eo = ei = e 2 ), which means that scattering 
by the boundaries does not take place, the Green function G s will be the Green 
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4-7T ||r — rg 



e iK e \\r-r \\ 



(40) 



where K 2 — e e Lu 2 /c 2 ac , and equation lp?T|l becomes the usual equation used in 
scattering theory by random media |171 1341 1351 I36| . In taking into account the 
boundaries, we have to change this Green function for an infinite random medium 
by Green functions taking into account the scattering by the boundaries. However, 
it is worth mentioning that the potential V docs not depend on the boundaries, 
but only on the random medium. Because of this property, we will apply exactly the 
same procedures developed in scattering theory by an infinite random medium, where 
the Green function for an unbounded medium must be replaced by Green functions 
describing scattering by boundaries. 

4. Link between the Green functions Gg 1 and scattering operator 

In this section, we show how to express the Green functions G^ 1 with the help 
of scattering operators which are common tools in scattering theory by rough 
surfaces These operators describe the field scattered by a rough surface 

illuminated by an incident plane wave. (The Green functions Gg describe the same 
phenomenon for a spherical incident wave.) For a rough surface separating two semi- 
infinite homogeneous media with permittivities eo and e e (Figure|3Jl , the wave equation 



z 



ey(Po) 




min x h{x) 



max x h(x) 



medium : £q 



e H (p) 



medium 1 : e e 



Figure 3. Incident plane wave with wave vector fc^ reflected into medium 1 
with the wave vector and transmitted into medium 1 with wave vector fcp~ . 
The polarization basis are depicted. 



can be simplified and transformed into the Helmholtz equation : 
(A + K$)E°(r) =0, z>h(x), 
(A + K*)E 1 (r)=0, z<h(x), 



(41) 
(42) 



with 
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and transversality equation : 

V-E°(x,z) =0, z>h(x), (44) 

y-E°(x,z) = 0, z</i(a;). (45) 

To find the fields E° and E 1 , we need the boundary conditions on the rough surface 
r s =x + h(x)e z : 

n s xE°(r s ) = n s xE 1 (r s ), (46) 

h s -e E°(r s ) =n s -e e E 1 (r s ), (47) 

n s x [V x E?°(r«)] = n s x [V x E 1 (r s )] , (48) 

n s • [V x £°(r s )] = n s ■ [V x E 1 ^)] , (49) 
and the radiation condition at infinity. For an incident plane wave 

E Ql (r) = E 0l (p Q )e lp °- x - iao(p o'> z 

coming from medium 0, the solution of these equations can be written on the following 
form [21]: 

E°(x,z) = E oi (x,z)+E 0s (x,z), (50) 
E 0s (x,z) = JAPR W (p\ Po ).E oi (p )e ip - x+iao ^ z for z>max x h(x) : 

(51) 

E l (x,z)=E u (x,z), (52) 
E lt (x,z)= J-f^T W (p\p )-E^( PQ )e^-^ a ^ for z < min x h{x) . 



(53) 



Here 



a (p) = ^ 2 -p 2 , a e (p)= VK?-P 2 . (55) 

It can be easily checked that the propagation equations (|41|) and l|42|) are satisfied 
with the representations I|51I53|) and the definitions l|55|l . To satisfy the transversality 
conditions P^)) and |@SJ, we need to decompose the scattering operators R (p\p ) 

and T (p|p ) on an orthogonal basis perpendicular to propagation vectors defined 
by p. These vectors are given by the following formula in medium and 1 (Figure 0J): 

k°+ =p + a (p)e z , (56) 

fcp~ = p - a e (p)e z . (57) 

|| In this definition, we need a precise the meaning of the square root because the integrand can 
be negative or complex if the media are absorbing. Since the imaginary part of the permittivity is 
always positive for an absorbing medium, we can use the following square root determination: 

^ z= (\A±^) 1,2 + i (\A^l) 1/2 , (54) 



which corresponds to classical square root operation for z £ '. 
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Figure 4. Wave vector decompositions 



The basis [e^~(p), ejy(p)] and [e v (p),e#(p)] respectively orthogonal to the 
vectors fcp + and fep~ are then dehned by 

- = e z X p , 



xk" 



o±i 



(p) 



\e z x 
e H (p) x fc" 
\e H (p ) x k ^ 



and 



eu(p) 



e z x k 



i± 



i±i 



(p) 



e#(p) x fcp* 



_+ Q o(P) 6 



IIPII 

#0 



IPI 



,\e H (p ) x fe^H i^e 
The scattering operators can be written with dyadic notations on these bases: 

R W (p\p Q ) = R w (p\Po)vv e° v + (p)e v (p )+R 10 (p\ Po ) HV e H (p)e y (p ) 

+i? 10 (p|p )y ff ey + (p)e ff (p ) + i? 10 (p|p ) H H e H (p)e H (p ), 

T 10 (p\p )=T m (p\p ) vv eir(p)e°-(p ) + T 10 (p|PoW e^(p)e°T(p ) 

+r 10 (p|p )y ff g£-(p)ejr(p„) + T 10 (p|PoW e H (p)e H (p ) , 

or in a matrix form: 

' # 10 (p|p Vv i? 10 (p|Po)^ ' 
i? 10 (p|p W ^ 10 (pIPo)hh 
T™(p\p ) vv T™(p\p ) VH 
T 10 (p\p ) HV T w {p\p ) HH 



R w (p\po) 



T W (P\P ) 



(58) 
(59) 

(60) 
(61) 

(62) 
(63) 

(64) 
(65) 



In a similar way, we can define scattering operators .R 01 and T" 01 which describe the 
reflected and transmitted fields when the source is in medium 1 (with permittivity 
e e ). For a rough surface z — h2(x) situated on the plane z = —H separating two 
homogenous media with permittivity e e and €2, we introduce scattering operators 
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R and T which describe the reflected field in medium 1 and the transmitted 
field in medium 2 when the source is in medium 1 (with the permittivity e e ). These 

21 21 

scattering operators can be obtained from the scattering operators R and T for 
a rough surface situated on the plane z = using the following properties of the 
scattering operators [2~T] : 

R H2 \p\Po) = e I (°"W+«.{Pb))ff i? 21 (p| Po ) , (66) 

T H21 (p\p Q ) = c^- a ^ +a ^ p "^ H T 21 (p\p ) . (67) 

In the rest of this paper, we suppose that we know the scattering operator expressions 
for R , T , R , T , R , T . Several approximate theories, like the small 
perturbation 37 , the Kirchhoff 20 , the small-slope approximation j^HI , the full- wave 
method [HHI , the integral-equation method 0H| , and others theories p3 |^ , 
can be used to obtain expressions for these operators |21|. With them, we can 
formally write the scattering operators for a slab with rough boundaries separating 
two homogeneous media (see Figure 0. We use the following notation S° (p| p ) for 

incident wave 




(2 



£0 






incident wave/ \ 


/ \ / l \ 
\ / S (p\Po) 







£2 



Figure 5. Scattering operator definitions. 

these scattering operators. The upperscripts a and b indicate the receiver location, the 
source location and if the waves are upgoing or downgoing. For example, the operator 
S 1+0 (p\ Pq) describes the amplitude of an incident downgoing wave from medium 
which is scattered into an upgoing wave in medium 1. The upwell electric field in 
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the medium 1 is given by 

E 1+ (x, *)=/t& ^■^^'S 1+ '°- (p\ Po ) ■ E 0i (p ) , (68) 

and 

(p\ Po) = S 1+Q -(p\p )w 4 + (p)4~(Po) + S 1+a -(p\p Q ) H v e H (p)e° v -(p ) 
+S 1+ °-(p\p Q )vH ev + (p)e ff (p ) + 5 1+0 -(p|p ) ff H e H {p)e H {p \m) 

_ , -r=ab . . rT -l(l —10 -=r01 —01 -=r21 —21 

io express the different operators b m function of it , T , it , T , it , T , 
we have just to add formally all the multiple scattering contributions on the rough 
boundaries ([2], p. 24). For example, S + (p\ p ) is given by 

■7=1+0- -=H21 —10 -^zH21 -=01 -=H2l —10 

S = R T +R R R T 

+R H21 -R 01 -R H21 -R 01 -R H21 -T m + ... (70) 
= R H21 -[I 1 - 1 --R 01 -R H21 ]-^T W , (71) 
where we have used the following notations: 

[A-B](p\Po) = /|^^(pIPi)-5(PiIPo)- (72) 

We have defined the projectors j^ 01 " by 

?i ""(P) = e 1 v a (p)e 1 v a "(p) + e H (p)e H {p) , (73) 
Ti ala °(p|Po) = (2^) 2 ^-p )7r a °(p), (74) 

where a and do are the sign + or — . The operator I ± °(p|Po) ^ s ^ ne unear identity 
mapping from the space vector defined by the basis [e.y °(p),eff(p)] to the space 
vector defined by [e.y (p),e.ff(p)]. To write the electric field E 1+ (x,z) inside the 
slab with the scattering operator S , we have implicitly assumed that the layer is 
sufficiently thick in order to have the condition — H + max x hi{x) < z < min x h\(x) . 
Furthermore, we will suppose that all the particles are inside the layer defined by 
—H + max x h2(x) < z < min x h\{x). Otherwise, the Rayleigh hypothesis |^ must 
be invoked to justify the use of the scattering operator for —H + maXxh^ix) > z and 
z > min x h±(x). 

Using the same reasoning for the different contributions, we obtain: 





= [I 1 - 1 --I? 1 .R H21 ]- 1 - 


li" 1 - . 


(75) 


s 1+1 ~ 


= R H21 .[f ± - 1 - + S 1 - 1 -], 




(76) 




= {l 1 ± - 1 - + S 1 - i -].R^\ 




(77) 


s 1+1+ 


= {1 1 + 1+ -R H21 .R° 1 }- 1 - 
= 7i H21 .S 1 - 1+ , 




(78) 
(79) 


s 0+0 - 


= 7i 10 + T 01 .S 1+1 -.T 10 , 




(80) 


s 1+0 - 






(81) 




= [I ± +S }-T , 




(82) 


s 0+1 - 


- T 01 ■ S 1+1 ~ 




(83) 


s 0+1+ 


= T 01 .[7i +1+ + 5 1+1+ ]. 




(84) 
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We now show how to express the Green operators G^ (r,r ) with the scattering 
operators S a a ° (p\ p ). 

To determine, for example, the Green operator G s (r,r ), we have to calculate 
the field produced inside the slab by a spherical source which is also in medium 1. We 
can decompose the function G s (r, ro) under the following form: 

S"(r, r ) = GHr, r ) + 1+ (r, r ) 

+Gl + 1_ (r, r ) + r ) + G*" 1+ (r, r ) . (85) 

The first term is the spherical source term in an infinite homogeneous medium (with 
permittivity e e ), and the following terms are the fields produced by multiple scattering 
process on the rough boundaries. We can obtain these contributions noticing that the 
source term 

„ /_ 1 \ giife ||J T-o|| 

Gr(r,r„)= J+-VV F.F.- - 

\ Kj J 47r||r — ro|| 

can be decomposed as a linear combination of plane wave using the Weyl formula 
jAppendix A| > : 

for z > zq : 



G^(r m) = - f d2p ° c iPn-(»-»o)+iq e (Pn)U-»o) (J 1 



(86) 



for z < z : 

G?(r, ro) = y^ e iPo-(*-* )- iQe ( Po ) (,-,o) (7 _ fe^^j _J_ ; (87) 

with fc p ^ = p ±a e (p ) e z and fc pQ = fc p ^/| |fc p ^ 1 1. Hence, for z ^ z , G^° is a linear 
combination of the following plane waves: 

which are transverse to the propagation directions fe po because (I — k po fe po ) • fe po = 
0. For each of these plane waves, we can calculate the field produced by the scattering 

1 X 1 X 

process on the boundaries with the scattering operators S . By using the 
superposition principle of the electric field 021 > we have for the contribution due to 
S 1+1 -: 

(J 1+1 ~( r rn \ — i f f d 2 P d2 Po ip-x-ip -x +ia e (p)z+ia B (p )z 

^^-2 J /i(2 7 r)2(2^) 2 

xS 1+1 ~(p\ Po )-(l-k po fe po ) 



(Po)' 

(89) 

- i - 

since in this case, the incident waves are propagating along the direction k p> . We 
summarize all the contributions l|85|l to G s under the form: 

Si 1 (r > r )=Sr°(r,ro)+ £ G^ la °(r,r ), 

(90) 
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with 

T< lala "( r ~\ — A f f_^P_^Pa_ ip-x-ip -xo+aia a (p)z-a ia e (p )z 

[r,r 0) 2 yy (27r)2(27r)2 e 

(p|p ) (J fepo fcpo . 

(91) 

where a, ao are the sign + or — . By using the same arguments, we demonstrate that: 
G s ° (r , r„) = G s + °~ (r, r ) + G s ~ °~ (r , r ) , (92) 
G°s(r, r ) = G° s + 1+ (r, r ) + G° s + l ~(r, r ) , (93) 
G° s °(r,r ) = G oo (r,r ) + G° + °"(r,r ), (94) 

with 

G la °~(r r n \ — - f [ d p d2p ° piP-as-ipo^o+oiaeCpJa+iaoCPoJ'o 
S [ ' ° j " 2 J J (2tt) 2 (2tt) 2 

x S la0 ~(p|p )-(7-fc Po fe° o ) 



"o(p ) 



(95) 



7; 0+la °|V ^ _ i /" / d?) d P pip-x-ip - a! o+iao(p) Z - go i Q!e (p ) Z o 

lr> r °)-2jj(2*)*(2*)* 

xS° + lao (p| p ) • (I - fep o ao fep °°) ' 



a e (p ) 



(96) 



=? 0+0_ r„ _ 1 f f d 2 Pp ip.x-ip -x +iQ (p)z+iao(p )zo 



<V (r,r ) 



(2tt) 2 (2tt) 2 



x5 0+0 -(p|p ).(T-CC^— 



"o(p ) 



(97) 



5. Lipmann-Schwinger equations and scattered field 



If we iterate equation l|31|l . we obtain the following series for the function G sv : 

G 1 S 1 V = G 1 S 1 + G 1 S 1 -V 11 -G 1 S 1 V , (98) 

= g* 1 + gI 1 -v 11 ■ gI 1 + <# • V 11 -cl 1 ■ V 11 ■ gI 1 + ... 

(99) 

If we introduce the transition operator T sv by 

Gsv = Gs + Z?s ■ T s l v ■ g" , (100) 

we obtain Lipmann-Schwinger equations in comparing the definition in (|100fl with the 
development in (|99|l : 

T 1 s 1 v = V n + V U -G s 1 .T 1 s 1 v , (101) 
Tlv = V U +Tlv -Cl 1 -V 11 . (102) 
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With these operators, we can straightforwdly rewrite in a compact form equations 
1271321 : 

(103) 
(104) 
(105) 
(106) 
(107) 

All the scattering processes in the random medium are contained in the transition 
operator T sv . If we know this operator, then we can calculate all the fields in the 
different media using (|103I107|I . Furthermore, for any source j source in the medium 
0, the incident electric field is 



7=700 

^sv 


7=700 




1 sv 


^s 


7=710 

^sv 


- ^s ' 




± sv 


^s 


^sv 


7=720 




1 SV 


—10 

^s 


7=f01 

<^sv 


7=701 


1- G°s 


± sv 




^sv 


= Gf- 




1 sv 


Gs 



E 0i (r) =iujn vac I d 3 r G (r,r )-j ( 



,(r ) 



(108) 



where [i vac is the vacuum permeability. The resulting field in medium is given from 
mm and by 

rrOs 



E 
E 



sv 

0s 
SV 



E 



0i 



E 



SV ) 



pOs , 7=r 



U S 1 SV 



'S I 



where 



£^ s (r) =iw/w / d 3 r Gs" (r,r )-j 



t=?°+o- 



£?5*(r) = iw/w d 3 r G s (r,r ) ■ j 



scarce 



rce( r o) 

(ro). 



(109) 
(110) 

(111) 
(112) 



The field -EgV is the scattered field produced by the random medium and the rough 
boundaries which is decomposed in two contributions where E s s (r) is the reflected 
field produced by an homogeneous slab (with permittivity e e ) with rough boundaries, 

and the second term G s ■ T sv ,7 ' 
and the rough surfaces. Vector E s l 
scattering by the particles. 



E s is the field scattered by the random medium 
is the transmitted field in medium 1 before any 



6. Ensemble average 

In order to calculate the coherent field reflected by the random slab, we need to 
define the averaging procedure. Let < • • • >s and < • ■ • >v denote, respectively, 
the ensemble average over the surfaces and the volume disorder. We also denote 
■C • • • ^$>sv the average over the rough surfaces and the volume disorder. We suppose 
that the rough surfaces and random medium properties are statistically independent, 
which means that <C • • • »sv=< [< • • ■ >s] >v=< [< • • • >y] >s- In the following 
development, we won't use any specific statistical properties of the rough surfaces 
(except that < h\{x) >s=< h 2 {x) >s= 0); thus, we don't need to specify them. An 
extensive description can be found in the references ^3 03] • The ensemble average 
over the random medium is defined by 

< f >v= d 3 ri...d s r N f(r u ...,r N )p(ri,...,r N ), (113) 
Jv x 
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where r 1 ,---,r^ are the particles positions, and p(ri, . . . , rjv) is the probability 
density function of finding the N particles at positions rx,---,rjsf- We will use a 
decomposition of this density function with conditional probabilities [161 117j : 

p(ri,...,r N ) =p{ri)p{ri,---,ri,---r N \ri), (114) 
p(n, ...,r N )= p{r i )p(r j \r l )p(r 1: ■ ■ ■ , f i: ■ ■ ■ f}, • • • , r N \r t , r 3 ) , (115) 

where the hat ""* indicates that the term is absent. The function p{r%) is the probability 
density function of finding a particle at n, p(ri\rj) is the conditional probability of 
finding a particle at Ti given a particle at rj, etc. If the particles are uniformly 
distributed inside the random medium V\, then the single particle density function is 
p{ri) = 1/Vi, where Vi is the volume of the area V\. In this case, we also define a 
pair-distribution function by 

9{\\r%-rj\\) =p{ri\r j )/p{r i ), (116) 

which depends only on the distance between the two particles if we suppose that 
the distribution of the particles is statistically homogeneous and isotropic. The 
normalization factor V\ is chosen in such a way that when the particles located at r% , 
Tj are far away from each over, their positions are uncorrelated (i.e., p(vi\Tj) = p(vi)), 
and we have 

Urn g(\\r i -r j \\) = l. (117) 

\\ri— t-j [ | — » + oo 

With these conditional probability functions, we define conditional averages: 



< / >v-.r l = I c\ i r l ...d?r i ...<Tr N f(r 1 ,...,r l ,...,r N ) 

xp(r t,..., fi,...r N \ri), (118) 
<f>v,r i ,r i = / d 3 r 1 . . . d 3 r, . . . d 3 rj ;. . . d 3 r N fin, . . . , r h . . . , v h . . . , r N ) (119) 

JVi 

xpin,.. .,fi,...,fj,.. .,r N \ri,rj) . (120) 



7. Coherent potential approximation and effective medium theory 

Until now, we have not clarified how to determine the effective permittivity e e (w). 
To determine e e (w), we use the fact that under some assumptions (mainly that the 
effective medium is not spatially dispersive fT7\ i.e., e(w,fc) = e(w)), it can be shown, 
using a diagrammatic technique, that the coherent part of the field < E >v which 
propagates inside an infinite random medium behaves as a wave in an homogeneous 
medium with a renormalized effective permittivity ^3 El H3 HH El 03 In order 
to insure this result in a self-consistent way, we introduce the Coherent Potential 
Approximation (CPA), which postulates that [TBIIT71 1271 15T1 OS) 

<G l sy>v=G l s. (121) 

This equation is in fact the definition of our effective permittivity e e and is a 
generalization of the classical (CPA) approach since we take into account the 
boundaries in the Green function definitions. Using equation \ 100(1 . we immediately 
see that condition 1|121J) is equivalent to 

<Tf v > v =0, (122) 
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where 

2& = V n + ^.<g.T£,. (123) 
Therefore, equations i|122M123|) provide a closed system of equations on the unknown 
e e (w) which takes place in the definition of V and G s . Because the random medium 
is made with spherical particles, it is convenient to express equation (|123|) as a function 
of the scattering operator t r for one particle located at r.; in a infinite medium 
| |Appcndix B| . It is defined by 

^^[i-vll-GTr'-vll (124) 
where the scattering potential U* 1 is: 

K] (r,r ) = (2tt) 2 S(r - r )< (r) , (125) 

<(r) = K 2 ac (e s -e 1 )Q s (r-r i )I. (126) 

In order to write equation 11231) with the operators t r ., following the Korringa 
demonstration |33| . we decompose 

F n (r,r ) = (2n) 2 S(r-r )K 2 ac (e v (r) -e e )7 (127) 

and 

T sv = V +V G s T sv (128) 
under the following form: 



where 

and 

with 



V U = V U + W U , (129) 



W (r,r ) = (27r) 2 S(r~r )K 2 vac (e e -e 1 )I, (130) 



Tsv — Tsv + Qsv i (131) 



Q S y = W 11 + W 11 • Gs ■ T S y . (132) 



By using the definitions (fT251 IT3U1 IH71 11251 we obtain 



A' 



i=l 

and in inserting (|129|) in i|123l) and using the definition l|131l I132|l , we have 
~ii ~ii ~ii — ii — ii 

T sv = V +V G s T sv (134) 
— ii —ii — ii —ii — ii 
= V + V G s (T sv ~ Qsv) • (135) 

In combining equation i|133|) and l|135|) . we decompose T S v m the following form: 

N 

Tsv = E C'sv,^ > (136) 

i=l 

where 

AT 

C" y ,„ = ^ + «ij • G" • (E ^,r, ~ Qsv) ■ (137) 
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If we subtract v* 1 ■ G s ■ C sv . r . in both sides of equation (|137fl . we obtain 

N 

(I-vil-G^-Csy.^^vll+vll-G'si V Cf v . rj -Q l s V ), (138) 



3=1, 



and then, 



where 



JV 

C SV,n = *5,r, + *5,r, ' ' ( CsVt, — Qsv) > (139) 

(140) 

= ^+^.G"-*L, (141) 

is the scattering operator for only one particle located at r, inside the slab. Using the 
decomposition in lj9"U|) : 

g" = GT+JG", (142) 



where 



rfe^ £ Gf". (143) 

a,ao=± 



In definition (|141(L we have 



[I v% ■ GT] ■ ts. r% = K- + Kl ■ SGs ■ ts >rt , (144) 



and we demonstrate that 
F 11 - f 



T ll T ll T ll r -^ n t11 /1 i f i 

tc Pl = t ri + t ri • 6G S ■ t s . t , (145) 



where 



^[7-^.GT]- 1 -^, (146) 



which is equivalent to 

F^^+^-GT-^. (147) 

The operator t* 1 describes the scattering process by a particle in an infinite 
homogenous medium (Figure EJ). We must be careful in determining the permittivity 
of this particle. In fact, as the propagator between two scattering events inside the 
particle is G 1 , the homogeneous medium surrounding the particle has the permittivity 
e e due to the definition of G^° | |Appendix B| . However, t r . doesn't describe the 
scattering by a particle of permittivity e s inside a medium of permittivity e e . If this 
was correct, the operator tjj, 1 defined by (|125l I12tjfl would contain the factor e s — e e 
and not the factor e s — t\. Thus, we have to renormalize the particle permittivity in 
introducing e s = e s — (ej — e e ), such as 

<A r ) = K U?' ~ ei)©s(r - njl, (148) 
= K 2 vac (e s -e e )Q s {r-r i jl. (149) 

Hence, the operator t r . is the scattering operator for a single particle of permittivity 
e s surrounded by an infinite medium of permittivity e e . 
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The operator t l g r . describes the scattering process by a particle located at 
inside the volume V\ of our slab (Figure EJ. If we iterate equation (|145fl . 

tl ri = t r ) + t r ) • 6Gf • t r ) + £ • ■ £ + (150) 

we see that the first term is the scattering process due to the particle, and the following 
terms describe the interaction between the particle and the boundaries (Figure|5J) since 
the terms SG S come from the slab surfaces. 

If we now look at equation l|139fl . we see that it describes multiple scattering 
process by different scatterers inside the slab. In fact, if we had defined the Green 
function G s without introducing the effective medium e e (w) but in taking the 
permittivity ei, the operators W and Q sv would have been zero operators (since 
e e (u>) = ei(w) in the definition (|130II132|) give null contributions), and the iteration of 
equation (|139[) show us the multiple scattering process inside the slab |271 I45| : 

N N N 

7=? n _ tH \ - -11 7=? n t 11 \ -> \ -ii —li-ii —li-ii . . 

^ SV,r t — t S,r i + Z-^i t S,r i 'Gs '*S,r 3 + ^ ^ t S>r . G S -tg^ G S 'tg^ +. . . . (151 J 

The operator C SVr . represents the field scattered by a particle located at 7"; 
which takes into account all the interaction effects with the other particles and the 
boundaries. In introducing the effective medium e e {uS) in the definition of G s , we 
see in equation (|139|l that the multiple scattering contributions are attenuated by the 
factor Qsv 

If we now average equation i|131|) . using the definitions 1113211130)1 and the (CPA) 
hypothesis (|122fl . we obtain 

(2tt) 2 S(r -r )e e K 2 vac I = (2tt) 2 S(r - r„) e x K 2 vac 1+ < f S y (r, r ) >v , (152) 

and since the ensemble average can be decomposed with the conditional probabilities 
I114|) and definition (|118fl . we have 

N 

< f]! v > v = < Cfv, ri >v , (153) 

i=l 

N r —ii. 

= V / d 3 rip(ri) < C sv > V;ri , (154) 
d 3 r, < C s 1 v - rj >y>, , (155) 

'Vi 

where we have introduced the particles density n = N/Vi- To obtain equation (|155|) . 

we have used the fact that < C sv . r . >y ;ri =< C sv . r . >v;r 3 f° r * h since we 
consider a statistical homogeneous random medium. 

In averaging equation (|139fl with the conditional average < . . . >y ;ri , and using 
definitions i|115l I120|l we obtain: 

N 

rii \ / ,q , -=n -=n 



j=l,j^i-^l 

< t Sr . ■ G s ■ Qsv >V;r, ■ (156) 
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Since . is the scattering operator for one particle located at rj , it only depends on 
the variable 7% and not on rj with j ^ i, and the average < • ■ ■ >v-.ri doesn't act on 
tg 1 ,,.. Furthermore, the averaging of equation l|132l) is 

<Q sv>v . r .= W + W ■ G s ■ < T sv > v , ri . (157) 

This equation is simplified by the (CPA) condition < T sv >v— 0, which can also be 
written under the following form: 

/ d 3 r iP (n) <r" y >y :ri =0. (158) 

As the identity in Ijl58|) is valid whatever the volume V\ is, we have < T sv >y- ri = 0. 
Thus, from equation (|157|l . we have: 

<Q 1 sv >v; ri = W U . (159) 

From the definition (|131I132I) and the coherent potential approximation, we deduce 
that: 

< Qlv >v = W n + W 11 ■ <?s ■ < Tf v > v , (160) 

= W U . (161) 

and 

< T l s V > v =< T S y > v + < Q~f v > v , (162) 

= < Qfv >v ■ (163) 
Accordingly, from the results H159U16llll63|) and equation (|155|) wc get for i 6 [1, JV]: 

< QSV >V;r t =< Tgy > V (164) 

= n / d 3 r, < Cly > V;rj , (165) 



and equation (|156|l is 

< Csv,r t >v-, ri = t s ^ + nt s G s ■ d J r 



J 



7=fH 



- < C S y rj >V; rj \ , (166) 

where we have used the approximation n ~ (JV— 1)/Vi which is valid for a large number 
of particle (N ^> 1). Using the same procedure, we can average equation l|139|) with 
the conditional average < . . . >y. ri r . and obtain an equation on < C SVr . >y ;ri ,r J 
in a function of < C 1 J Vrk >y. ri r . rk and so on. Hence, a hierarchical system of 
equations can be generated on the unknown < C SVr . >y ;r; , < C SVr . >v-.r i ,r i - I 

< C ' sv.r k >v-,ri.rj.r k ,- ■ ■ ■ We close this system by using the Quasi-Crystalline 
Approximation (QCA) which states that [Tfil IP71 l8H \M HH] 

• C SV,r \ r ,r C SV,r 'V;r ■ (167) 

As was demonstrated by Lax this approximation is strictly valid when the 

particles have a fixed position, as in a crystal. The quasi-crystalline approximation is 
equivalent to neglect of the fluctuation of the effective field, acting on a particle located 
at rj, due to a deviation of a particle located at from its average position. Under 
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this approximation, the effective permittivity e e (aj) satisfies the following system of 
equations: 



(2%y5(r - r ) e e K l vac I = (27r)^(r - r ) d K ac I + n d d r< < C SKr . (r, r ) >y :ri , 

(168) 



^11 



< ^SKr; >V>\ — t S r . + 71%^ ■ Org • 



dSkGki-i-iiD-i] <c 



Vi 



SV.r, 



>V;r 



(169) 



We can simplify these equations by noticing that the contribution 8G S due to the 

boundary in G s = + 5G S can be neglected in equations (|169fl and i|145[) if the 
following condition K'JH ^> 1 with K'J — ImK e is satisfied. Usually, we define the 
extinction length as l e = l/2K'J and we see that the previous condition means that 
the slab thickness must be greater than the extinction length. 

For example, if we analyze the contribution G s of 5G S , we have 



77I+0- 



-A [[ &2 P0 Jp--iPo-«o+i 



.(p) 2+ iao(Po) 2 o 



(2tt) 2 (2tt) 



(p| p ) • (J-fc° fc° ) 



"o(Po) 



(170) 



where 

S 1+ °~(pIPo) = 



d 2 Pi d 2 p 2 + i( ae (p) +ae ( Pl )) g 
(2tt) 2 (2tt) 2 



X R (p\pi) 



-=1 — 1— , . , -7=1-1- 

*x (Pi|P 2 ) + S 



(P1IP2) -r'ValPo) 



where we have used the following property: 



R H 21 (p |pj) = e +i(«e(p)+ Qc (p 1 )) # 21 (p | pi ) 
We see that ° contains a factor: 

g ia e (p) (z+-H")+iao(p )zo 

which is negligible far from the lower boundary if we have K"H ^> 1, since 



0. 



(171) 
(172) 

(173) 
(174) 



Similarly, we show that far from the boundaries, the other contributions to 5G S are 
negligible compared to G l . Thus, we replace in equation i|169fl the term G s by G 1 
and the operator t s r . by t r . , and we have 



< C 



SV,r, 



11 



vll 7=?°° 
nt r . • G : 



d 3 rM\\rj ~ n||) - 1] < C sv , Tj > V;rj (175) 



In doing this, we neglect all the boundary effects in the calculation of the effective 
permittivity e e , and equations obtained are the same used to calculate the effective 
permittivity in an infinite random medium |16U17ll55] . In an infinite medium, we can 
use the Fourier transform to write the equation l|175fl , which is defined by 



/(fe|feo) = J d 3 r d 3 r exp(-ifc • r + ik Q • r ) f{r, r ) . 



(176) 
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In the Fourier space, the translational invariance of the infinite medium implied the 
following property of the scattering operators t r . ( |Appcndix B| ): 



t (fe|fc ) = exp(-i(fc - fe ) ■ n)t (fe|fe ) . 



(177) 



where t]^(k\ko) = t r 1 =o(^ ; l^ s o) is the scattering operator for a particle located at the 
origin of the coordinate (see 1171 US) ). Using the property of (|177|) and equation 
l|175|) . we show that < C sv r . (fe|fc ) >v-,n verifies also a property similar to 11771) : 

<C 1 s 1 v ^ ri (k\k ) >y ;ri =exp(-i(fc-fc )-r. i )c"(fc|fe ), (178) 

>y ; r i= o- By using the 



< C«T/, n(fc|fco) 



where we have defined C a (fe|fco) — - ^svr i= o 
properties (|177U178l) in equations i|175ll68fl . we obtain 

e e K 2 vac I = a K 2 vac I + n c" (k | fc„) , 



C Q (k\k ) 
where 



t a (k\k Q )+n 



) , (179) 
- fc x ) (^(fclfci) ■ Gr(fei) ■ ^"(fci|fco) ,(180) 



—11 11 11 -3=700 —11 

v^ir^o) = S(r - r )tT"(r) , 
v 1 1 (r)=K 2 vac (i s -e e )Q s (r)I 



(181) 
(182) 
(183) 



and 



h(k-kx) = J d 3 r exp(-i(fc-fci)-r)[g(||r||)-l], (184) 

GT(fe) = y d 3 r exp(-ife • r)G^(r) . (185) 
In equation l|185|) . we have used the translational invariance of the Green function: 




G™ (r,r ) = G 1 (r - r ). Formula (tT79l TO is a non-linear system of equations 
on the unknown e e (oj). If we neglect the correlation between the particles (i.e., 
h(k — k\) = 0) and define the Green function G 1 in replacing the effective permittivity 
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e e by ei, we obtain the Foldy's approximation also called the independent scattering 
approximation (ISA) [271 E7J H3 HI El : 

e e K^ ac I = e 1 K^ ac I + nt o 1 (k \k ) • (186) 

However, this result is greatly improved under the (CPA-QCA) approach since for 
Rayleigh scatterers, an approximate formula for e e (w) can be derived from equations 
<|179l I180|) , which is a generalization of the usual Maxwell-Garnett formula El 1481 
1491 15U) . One can also obtain an approximate formula for the effective permittivity 
which at the same time contained the Maxwell-Garnett formula and the Keller 
approximation |51j . 

8. Coherent field 

By using the expression in I|11U|) and the (CPA) condition (|122|) . the average electric 
field is 

< E° sv > v = E° l + E° s s . (187) 
If we average over the surface disorder, we have 

« E° sv » sv = E M + < Ef > s . (188) 
Hence, for an incident plane wave, 

E 0l (x, z) = E 0l (p ) e'Po^-^olPo) ^ ( ( 189 ) 

the field in the medium is 

« E° sv (x,z) » S y= Ef*(p Q )J*>"- ia °<*>)', (190) 

+ / -^j p - x+iao{v)z <S° + °~(p\p )>s -E^po), 
where equations fTBI I7HI l%U)l give us 

S° + °- =R W +J fil -R H21 ■ [I 1 - 1 --R Q1 -R H21 ]- 1 T W . (192) 
For statistical homogeneous rough surfaces, we have [2"T] : 

< S° + °-(p\p ) > s = (27r) 2 S(p-p )S Coh (p a ), (193) 

where S (p ) is a diagonal operator: 

S c ° h (p Q ) = S c ° h (p ) vv e Q v + ( Po )e ( l-(p )+S coh (p Q ) HH e H (p Q )e H (p ) , (194) 
or in a matrix form: 

p-M - ( s *°y™ ^ ) , < 195 , 

and 

< E° sv (x,z) » w = E m (p )e ip °- x - iao{p ° )z + S coh (p a ) ■ E°\p a )e ip «- x+ia "^ z . 

(196) 

Hence, the coherent field behaves as if the slab was an homogeneous medium of 
permittivity e e with planar boundaries but with modified Fresnel coefficients given 
by the two diagonal elements of S^ at (p ). 



(191) 
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9. Application 

Few approximate theories give explicit expression for the scattering operators H192[l 
for a slab. Most of them use the small-perturbation theory [S3 E3 \EM EH1 EE ED 
to derive some approximate expression of S , but the perturbative development 
needs to go up to the second order to take into account the roughness of the surfaces. 
We must also mention that the Kirchhoff theory in the geometrical optics limit and 
the full- wave method have been extended for a slab with rough boundaries [SSJ 1513 • 
However, we know that for the coherent part of the scattered field, the exponential 
term present in the Kirchhoff theory gives an accurate description even in the small- 
perturbation limit IS!]]. Thus, under the Kirchhoff theory (or the first term of Small- 
Slope Approximation, which gives the same results EU), we have for Gaussian rough 

surfaces eh nana, 



< R 10 (p\ Po ) >s= (2n)H(p-p )r w (p )e- 2a ^ 2 ^ 



< T u (p|p ) > s = (2tt) 2 5(p - p ) r (p )e- 



-("o(P )-"e(Po)) 2 ^?/2 



-H21 



< R (p| Po ) > s = (2nY 8{p - p ) r H 21 (p )e 



-01 



< R (p|p ) > s = (2tt) 2 5(p-p )r '(p )e 

< T 10 (p\p ) > 5 = (2n) 2 S(p-p )t W (p )e- 



where <j\, o~i are the rms- heights of the rough surfaces: 



-2a e (p ) 2 <xJ 

5 

-("e(p )-"o(Po)) 2 ^?/2 



(197) 
(198) 
(199) 
(200) 
(201) 
(202) 



and r ,t ,r 



0"! 



H21 



V< h x {x) 2 >< 



0-2 = \/< h 2 {x) 2 > s ■ 



(203) 



r 01 , r 01 are reflection operators for the planar surface. We write 



these operators with matrices. For example, 

r W (p )=r w (p ) vv e Q +( P() )e v -(p )+r w (p ) HH e H (p )e H (p )(204) 



is written 



and we have: 



[r w (Po)] = 



,.io 



(Po) 



vv 







r w (p Q ) HH 

e eno(Po)-eoae(p ) 



^ iU (Po): 

*%o) 



e e Qo(p )+coae(p ) 





a a(Po)- a e{Po) 
a o(Po)+ a <=(Po) 



2(e e e )2 q e ( Po ) 

; eQ!o(Po) + «Oae(Po) 





2a e (Po) 

Q o(Po) + "e(p ) 



\r H21 



\r 01 



(Po)} = 
(Po)]=-[r W (Po)}, 



£2ae(p )-e e «2(p ) 

£2Qc(Po) + eeQ2(p ) 





«e(Po)-«2(Po) 
a e (p ) + a 2 (p ) 



3 2ic* e (p )H 



* 10 (Po) 



2(eoe e )? ao(Po) 

C e Q o(Po) + fO«e(p ) 







2 "o(Po) 
Qo(p )+ae(p ) 



(205) 

(206) 

(207) 

,(208) 
(209) 
(210) 
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Hence, in using the independent scattering approximation for the rough surfaces to 
calculate < S , ° + ° >s, 

< S° +a - > s =< R 10 > s + < T 01 > s • < R H21 > s (211) 
■[t ± - 1 --<R 01 >s-<R H21 >sr 1 -<T W > s , (212) 
we obtain an approximate expression for the diagonal matrix [S ,c °' i ] given by 
[S coh ( Po )} = [r w (p n )}e- 2a «^ 2 ^ +c -(«a(P )- ac (P )) 2 *l-2^( Po f*l pOi (po)] _ [¥ H2i {po)] 

. ([j 2 ] _ e -^e(Po) 2 K + ^) [F oi (po)] . [r^ 21 (p )]) _1 • [t 10 ( Po )] , (213) 

where [J 2 ] — diag(l, 1) is the two dimensional identity matrix. Using the following 
identity [r 10 (p )] 2 + \i (p )] ■ [t (p )] = [1 2] ( which is the conservation energy law 
for a planar surface) and the property in 1)209(1 . we rewrite equation (|213|l under the 
following form: 



[S (po)] = [[c\p )] ■ [r 10 (p )] + [c 2 (p )] • [r 11 21 (p ) 

•([7 2 ] + [c 3 (p )] • [r 01 (p )] ■ [r H 21 (p )]) * , (214) 
with 

[cHPo)]=e- 2a ° (p ° )2ff? [7 2 ], (215) 
[c 2 (p )] = e _ ( Qo ( po ) _Q ^ Po ^ 2 °"?~ 2a <=(Po) 2 CT 2 f[j 2 ] — ( e ~( Q °( p o) +ac ( p o)) 2 CT i - l)[r 10 (p )] 2 ^ 

(216) 

[c 3 (p )] = e ~ 2tt =( p o) 2 ( <7 ?+ <7 =) [7 2 ] . (217) 
For a random medium with planar boundaries (ctj = <t 2 = 0), we obtain 

[S^CPo)] = ([r 10 (p )] + [r ff21 (p )]).([7 2 ] + [r 10 (p )].[r ff21 (p )]) _1 ,(218) 

which is a diagonal matrix which contains the usual reflection coefficients for a planar 
slab separating three homogeneous medium with the permittivities eo,e e ,e 2 . (See 
references ^3 El E21 I53|.) In comparing expression (|214|l and (|218|l . we see that 
the rough surfaces modify the reflection coefficients for a planar slab by adding new 
factors [c ],[c ],[c ]. The random medium doesn't change the form of the reflection 
coefficients but only the permittivity t\ of the initial medium by an effective one e e . 
Furthermore, if the random medium is highly scattering and thick, the imaginary part 
of the effective permittivity e e is important, and the factor exp(2 i a e (p ) H) in the 
expression (I2()8|l of [r^ 21 ] is very small. The contribution [r^ 21 ] in equation l(218|l 
becomes negligible compared to [r 10 (p )], and we have 

[S coh (p )} = e- 2a ^o) 2 ^ [r w (p Q )} , (219) 

which is the Kirchhoff term for a slab separating two semi-infinite media with the 
permittivity eo and e e . In this case, the lower boundary doesn't contribute to the 
coherent field. 
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10. Conclusion 



We have considered the scattering of an electromagnetic wave by a random medium 
with rough boundaries. We have formulated the solution of this problem using two 
kinds of Green functions. The first one describes the scattering by the rough surfaces 
and the random medium, and the other represents the scattering by an homogeneous 
slab with rough boundaries. As equations obtained are similar to those used in 
scattering theory by an infinite random medium, we were able to introduce the 
coherent potential with the quasi-crystalline approximation to calculate the effect of 
the random medium on the coherent field. With this approach, the random medium 
contribution is taken into account by an effective medium permittivity. The surface 
scattering contributions on the coherent field are included in the scattering operator 
of the system, which describes the scattering by the rough boundaries. This operator 
can be approximated using the usual scattering theories by rough surface like the 
small-perturbation, the Kirchhoff, or other more sophisticated theories. To derive 
these results, we have supposed that the slab is sufficiently thick to insure, for one 
hand, that their exist a layer (-H + maxxh-iix) < z < min x hi(x)) between the two 
rough boundaries which contains the scatterers, and in second hand, that the effective 
permittivity e e (ui) doesn't not depend on the boundaries (K"H ^ 1). 

In the following papers, we will use our Green function formulation of the 
scattering problem to derive a radiative transfer equation describing the scattered 
incoherent intensity. Furthermore, we will investigate the case of an highly scattering 
medium where a vectorial diffusion approximation permits simplifying the radiative 
transfer equation. 



Appendix A. Green functions 

Appendix A. 1. Scalar Green function 
The solution of 

(A + K%)Go(r,r ) = -5(r-r Q ) (A.l) 

in an infinite medium which satisfies the radiation condition at infinity is a generalized 
function given by 

piK ||t — r \\ 

where P.V. is the principal value defined by 

P.V. [ d 3 r Gg°(r - r ) <f>(r ) = lim / d 3 r Gg°(r - r ) 0(r o ) , (A.3) 

J a ^°JV-V a (r) 

where <p(r) is a test function, V a (ro) is an exclusion volume with size a around the 
singularity located at r. In equation l|A.2|) . the exclusion volume is a sphere [ST]. This 
generalized function can be represented as the usual spherical function for r/rj, 

pi K | |i — t-o 1 1 

G^(r-r Q )= (A.4) 
47T \\r - r || 

Using Fourier transform and the residue theorem, we can also write this function 
under the following form: 

1 f Jp-x+iao(p) \z\ _ 



G ° {r) ~ 2 J (2tt) 2 6 a (p) r ^°' (A ' 5j 
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The general solution of equation (|A.1|I can be expressed with equation (|A.5|I as a 
generalized function: 



G^{r)=P.V. 1 - f-!pL c ip-*+ia„( P )W 1 r ^ Q 
2 J (2tt) z a {p) 



(A.6) 



where the exclusion volume used in the definition of the principal value (|A.3(1 is a 
pillbox of arbitrary cross section but thin in the z direction due to the term \z\ in the 
expression l|A~6|l of the Green function [52llH5| . 

Appendix A. 2. Dyadic Green function 
The solution of 

V x V x G °>,r ) - Kl G °°(r,ro) = S(r - r )7 (A.7) 
in an infinite medium is a generalized function given by: 

G °> - r„) = (l+ ^VV) G °°(r -r„) „ (A.8) 



which is short notation for 



(A.9) 



J d 3 r G °°(r - r ) <f>(r ) = (i + ^wj J d 3 r G~(r - r ) 0(r o ) , 

where VV/(r) = V[V/(r)]. When we use the representation (jA. 2|> or l|A.6j) in (|A.8j) . 
the action of VV on the exclusion volume V a (r) produces a singularity: 

G °°(r - r„) - P-V. (i + ^Vv) G^r - r ) - -^<5(r - r )X, (A.10) 

where the principal value is defined by 



P.V. J d 3 r 0(r o ) f I + ^ vv ) G S> - r ) 



= lim/ d 3 r o 0(r o ) f 7 + Iw) G»(r - r ) . (A.ll) 

The operator X depends on the exclusion volume chosen; for a spherical volume, we 
have L = J/3, and for a pillbox thin in the z direction, we have L = e z e z . As r ^ Tq 
in the principal value term, we can use the representation (|A.5|I to calculate the first 
term in l|A.10fl and we obtain 

Gn°°fr m) = -PV ( d2PO C iPn-r»-»oHiao(Pn)U-ZQl (J /; 0± /; 0± -) I 



1 



-^(r-r )J, (A.12) 

-0± 



where k p = p + sgn(z — Zq) ao(p ) e z and 

+ 1 , if z > 

sgn(z) = { if 2 = . (A.13) 
-1 if z < 

The upperscript sign in k is given by the sign of the function sgn(z — zq). 
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Appendix B. Transition operator for one scatterer 

The electric field produced by an incident wave E l (r) scattered by a spherical particle 
of radius r s , located at rj, with a permittivity e a , and surrounded by an infinite 
medium of permittivity tb is given by 1171 146j : 

E{r) = E\r) + J d 3 n T2?{t,t x ) -^.(n) • E{r{) , (B.l) 



where 



and 



vL (n) = K 2 vac (e a - e b ) 6 s (n -r,)I, (B.2) 



6s(r) -^ if ||r||>r. ' (R3) 



The Green function G b (t*,t*i) is defined by 



Gf - r ) = (l+ ^VVj Gg°(r r ) , (B.4) 

gi-R"b II' — roll 

G?(r-r ) =P.V , (B.5) 

47T | |r* — t*o 1 1 

where K% = ebK% ac . The transition operator for one particle is defined by: 

E(r) = E\r) + J d 3 n d 3 r 2 G™(r, n) • i"(n, r 2 ) • ^(r 2 ) . (B.6) 

In comparing the definition in (|B.6|) with equation (|B.lf) . we obtain 

t"(r,r ) = ^.(r)<5(r-r )+ J d 3 n ^(rJ.GTCr.nJ-t^n.roJ.CB.T) 

In the Fourier space, we have 

t r 1 .(k\k ) = v 1 r .(k-ko)+l |^ ^.(fe-fe 1 )-GT(fe 1 )-t^(fe 1 |feo),(B.8) 



with 



and 



(2tt) 3 

vl. (k-k ) = J d 3 r e - i ( fe - fe °)- r v x r . (r) , (B.9) 



Gr(fei)= J d"rc-^ r G,(r). (B.10) 
We easily check that tjj, . verifies the following property: 

vl 3 (k-k ) = e- 1 ( fe - fe °^ vl(k - ko) , (B.ll) 

where Uo(fe-feo) = v]. =0 (k — ko). In iterating equation (|B.8(1 . and using the property 
(|B. 1 1|) . we demonstrate that 

t"(fe|fe ) = exp(-i(fe - ko) ■ r 3 )t o \k\k ) , (B.12) 

where i* 1 = t r . =0 is the transition operator for a particle located at the origin of the 
coordinate. If we consider an incident plane wave E l (r) 1 

E l (r) = E t (k )e ik °- r , (B.13) 
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transverse to the propagation direction fc : 

k o -E l (k o )=0 ^ (l-k ko)-E l (k ) = E l (k ) (B.14) 

where fco = Kb fco, and fco • fco = 1. The far-field scattered by a particle located at the 
origin is given by 

# 6 ||r|| » 1 , E(r) = E i (r)+ ^ Kb !. M ! (7-fcfc)^ 1 (k\k y(l~k Q k yE l (k ) , (B.15) 

47T ||r|| 

where 

t"(fc|fc ) = J d 3 rd 3 r e- ik r+ik " r "t 1 o 1 (r\r ). (B.16) 
We have used the following far-field approximation to derive the equation (|B.15(1 : 

G b (r,r 1 )«^ w (J-fcfc)e- lfc - 1 , (B.17) 
47r ||r|| 

where k — K\,r. Usually, the far field scattered by a particle is written in the following 
form g7l El ESJ: 

Jf 6 ||r||>l, E s (r) ^ ^^y 1 J(k\k ) ■ E l (k ) . (B.18) 

IMI 

For spherical scatterer, an exact expression of this operator is well-known and given 



by the Mie theory 47, 64, 66, 65 . In comparing equations (jB . 1 B|> and (|B.18J) . we have 
the following relationship between /(fe|fco) an d t a (k\ko) : 

47r7(fc|fe ) = {I - kk) ■ il\k\k ) • (7 - fcofco) , (B.19) 

where k — K^k, and fco = -Kfcfco- The operator t^ 1 (/c | fco ) is a generalization of the 
scattering amplitude /(fc|fco) since it contains also the near-field component scattered 
by the particle. Furthermore, we see that the far-field component is obtained in 
taking only the transversal components of t Q (fc|fco) (due to the projectors I — kk and 
J — fcofco) and using an on-shell approximation (since fc • fc = fco • fco = Kb). 

Appendix C. Reciprocity of C G (fe|feo) 

In the section [7] we have used the following decomposition of the operator Tgy: 

T s \ = V^ + V^.g\}-T s \. (CI) 
But we could have used this equivalent formulation: 

T^^+T^.G^.V". (C.2) 

. _ ~n 

According to the section Q we define a new operator V such that: 

V ll = V U + W l \ (C.3) 

where W U is defined by equation l|13(J|) . The operator T sv can also be decomposed 

~n n . . 

in defining new operators T sv or TT SV and using either equation (|C.1(I or (|C.2ll : 

with the Eq. m I Si = Si + Si ' 7,n =11 (C4) 
I Qsv = W +W G s -T sv , 

with the Eq. Q { TIM = ^i + ' ii W n (C5) 

{ QQsv =W +T SV -G S W , 
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Following the demonstration of sectional we derive the following equations: 

± SV — 2^,1=1 ° SV,r z ) trt „s 

r 11 -7 11 + 7 11 .r 11 .^ r 11 o 11 ^ ( ' 

^ SV,ri — z S,r z ' z S, ri <^S ^Z^ J = lj^i ° SV, rj ~ W SV 



- 11 



T SV - Si=l CCgV , 

rr 11 -7 11 +/v JV rr 11 oo'M.r 11 ./ 11 ( j 



If we define the average of the operators C SVr . and CC SVr . at the origin (r, = 0) 
by: 

C Q =< C S v ir . =0 >y ; r,=o , (C.8) 



CC a =< CC SVr , i=0 >v,n=o , (C.9) 

(CIO) 

we obtain, by using equations l)C.6IC.7|) . the following expressions: 

d 3 fci 



/A o 1 
— -| - fci)i"(fc|fcO -ST(fci) • C"(fei|fe ) 



(C.ll) 

Cc"(fe|fco)-*o 1 (fe|fc ) + n y h(k 1 -k )CCf(k\k 1 )-~&?(k 1 )-tt 1 (k 1 \k ), 

(C.12) 

Furthermore, under the (CPA) approximation we have < T sv >v = 0, and from 
equations (|C4IC5|) . we deduce that < Q sv >v=< QQsv >v= W and then 

< Tg V > v =< TT S \ > v . (C.13) 

By using the decomposition ljC.6l IC.7H and the definition of the conditional average 
< >v-n, equation (|C.13|) can be written as: 



/ d 3 r 4 <cf Vri > v , ri =n / d 3 ^ <CC s \ r . > v . rt 

JVi JVx 



(C.14) 



This identity is valid whatever the volume V\ and the position Ti of the scatterer, and 
thus we have: 



C a (k\k ) = CC o (fe|fe ). (C15) 
Hence, the operator c] ) 1 (k\ko) satisfies the following equations: 

c"(fc|feo)=t"(fc|fe )+n J T^yjj h(k - fci) t"(fc|fci) • G^°(ki) ■ (fci |fc ) , 
Cl\k\k ) = t"(fc|fc )+n J !^h{k 1 ~k )cl\k\k 1 )-G™{k 1 )-t o \k l \k ). 

However, since the operator 7* 1 is reciprocal and G^°(ki) = G^°(—ki), we easily 
show using equations (jC.lllC.12|) that C Q (fc|fc ) = [CC D (— ko\ — k)] T where T is 
the transpose of the operator. From the identity (|C.15J| . we conclude that the operator 

C Q (fe|feo) is reciprocal: 

C 1 \k\k ) = {C 1 \~ko\-k)] T . (C.16) 
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